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Abstract 

A multidimensional computational 
model of the pressurization process In a 
slush hydrogen propellant storage tank was 
developed and Its accuracy evaluated by 
comparison to experimental data measured 
for a 5 ft diameter spherical tank. The fluid 
mechanic, thermodynamic, and heat transfer 
processes within the ullage are represented 
by a finite-volume model. The heat and ma s s 
fluxes at the ullage boundary were computed 
In auxiliary analyses and specified as Input 
to the finite-volume model. The model was 
shown to be In reasonable agreement with 
the experiment data. A parameter study was 
undertaken to examine the dependence of 
the pressurization process on Initial ullage 
temperature distribution and pressurant 
mass flow rate. It Is shown that for a given 
heat flux rate at the ullage boundary, the 
pressurization process Is nearly Independent 
of initial temperature distribution. The mass 
flow rate study revealed decreasing pressur- 
ant mass requirement with Increasing 
pressurant mass flow rate. Further, signifi- 
cant differences were identified between the 
ullage temperature and velocity fields pre- 
dicted for pressurization of slush and those 
predicted for pressurization of liquid hydro- 
gen A simplified model of the pressurization 
process was constructed In search of a 
dimensionless characterization of the 
pressurization process. It Is shown that the 
relationship derived from this simplified 
model collapses all of the pressure history 
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data generated during this study Into a 
single curve. 


Nomenclature 

Cv Specific heat at constant volume 

hg Specific enthalpy of pressurant gas 

m g Pressurant mass flow rate 

p Tank pressure 

Q l/s Heat transfer rate to Uquid/slush H 2 
Q w Rate of heat transfer to the tank wall 
R Gas constant 

t Time 

t p Total pressurization time 

T Temperature 

V Tank ullage volume 

p Density of the gas In the tank 

Subscripts 
f Final state 

1 Initial state 


Introduction 

Vehicles which use cryogenic propel- 
lants because of their high specific Impulse, 
Incur the difficulties associated with manag- 
ing the propellant prior to its use. The 
propellant management issue of Interest In 
this study is the pressurization of a vehicle 
propellant storage tank prior to engine Igni- 
tion. This pressurization may be required to 
alleviate feed-line cavitation problems or may 
provide the sole motive force for transporting 
propellant from the storage tank to the 
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engine. Since all methods proposed for tank 
pressurization require on-board resources 
(such as pressurant gas), optimization of the 
pressurization process Is Important for effi- 
cient vehicle design and operation. 
Optimization requires an understanding of 
the complex fluid dynamic and thermody- 
namic processes occurring In cryogenic 
propellant management systems and 
research into these processes has been 
pursued for over 20 years 1,2,3 . 

The proposal to use a mixture of solid 
and liquid hydrogen (known as slush hydro- 
gen) as the storage medium for the National 
Aerospace Plane (NASP) presents a 
formidable challenge to understanding the 
pressurization process in the propellant stor- 
age tanks for such a vehicle. Modeling the 
pressurization process requires simultane- 
ous solution of fluid mechanic, thermody- 
namic, and heat transfer models for a tank In 
which three phases of matter are present. 
Most existing models 1,4 are fundamentally 
one-dimensional and based on simplifying 
assumptions and experimentally derived 
correlations. Although these models have 
served well In the past, their range of validity 
Is limited by the range of cases from which 
they are derived. It is not clear that they are 
adequate as design tools for future vehicle 
development. Since the flow field and tem- 
perature field In a tank during pressurization 
are in actuality multi-dimenslonal, it is rea- 
sonable to expect that accurate modeling of 
the pressurization process will require inclu- 
sion of the multl-dimensional effects. One 
study 5 , using the same computational 
technology as the present study, has been 
performed to examine certain multidimen- 
sional effects. The current study was under- 
taken to identify or develop a tool for multi- 
dimensional modeling of the pressurization 
process and to use that tool to begin a study 
of the parametric dependence of the 
pressurization process. 

The first phase of this study identified 
the FLOW-3D code 6 as a good candidate for 
modeling the pressurization process and 
demonstrated its accuracy by comparing 
computational predictions to experimental 
data. The second phase was a study of the 
parametric dependence of the pressurization 


process on initial ullage temperature distri- 
bution and pressurant mass flow rate . The 
third phase was development of a dimension- 
less characterization of pressurization by 
construction of a simplified model and exam- 
ination of the physical processes which occur 
In the ullage during pressurization. The sec- 
tions of this paper are organized to reflect 
this sequence of efforts. The Anal section 
summarizes the results of the entire study 
and presents conclusions drawn from these 
results. 


Experiment Configuration 

The experiment configuration modeled 
during this study was a 5 ft diameter spheri- 
cal tank partially filled with liquid/slush 
hydrogen. The tank was pressurized by 
Injecting gaseous hydrogen through a 1 ft 
diameter hemispherical diffuser located at 
the top of the tank. The pressurization pro- 
cess was Initiated at the appropriate satura- 
tion pressure, (1.1 psia for slush hydrogen or 
17.4 psia for liquid hydrogen), and was 
terminated when the tank Internal pressure 
reached 50 psi. Although the experimental 
procedure then proceeded to a hold period 
followed by expulsion of the propellant, the 
focus of this study Is on the pressurization 
process. A schematic of the tank configura- 
tion Is shown In Fig. 1. A typical tank pres- 
sure history recorded during the pressuriza- 
tion of liquid hydrogen Is shown In Fig. 2. 


Evaluation of The Computational Model 

The first phase of the study required 
Identification or development of a computa- 
tional tool capable of modeling the pressur- 
ization process. At a minimum, this requires 
the ability to model the transient low-speed 
compressible flow of an ideal gas, Including 
heat transfer effects, in a reasonable repre- 
sentation of the geometry of interest. An 
important secondary consideration was the 
ability of the tool to also model the expulsion 
process. Since the fluid level in the tank 
changes during the expulsion process, this 
requires modeling the flow of an Incompress- 
ible fluid with a vapor/liquid Interface that 
changes position with time. The FLOW-3D 
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code 6 was selected as the highest rated 
candidate because it appeared to meet all of 
the basic requirements and because of previ- 
ous successful applications to similar prob- 
lems 5 - 8 . 

FLOW-3D Is designed to model tran- 
sient flows of fluid and thermal energy In 
three space dimensions In complicated 
geometries which may Include free surfaces. 
Either Cartesian or cylindrical coordinates 
may be used to define the problem geometry 
and two-fluid calculations can be performed. 
Some of the models and features which con- 
tributed to selecting FLOW-3D for this study 
include: a compressible flow model (Ideal 
gas), a two-fluid interface model, time- 
dependent boundary conditions, a thermal 
buoyancy model, a surface tension model, 
and a turbulent flow model. To efficiently 
use the available computational resources, 
the code is divided Into a pre-processor, a 
main-processor, and a post-processor. The 
pre-processor resides on a VAX and Is the 
interface through which the problem was 
defined. The data is then transferred to a 
Cray XMP where the main processor 
performs the analysis and generates the 
output flies which are returned to the VAX 
for post-processing. A typical run for this 
study, using a two-dimensional pie-sector 
model of the tank, required approximately 
25 minutes of CPU time on the Cray. 

The first step In the evaluation process 
was to determine a suitable computational 
model for the spherical tank. Exper imental 
data 1 obtained for pressurization of a tank 
with 56% ullage by volume was selected to 
test the performance of FLOW-3D. The free 
surface of the liquid hydrogen was modeled 
as a solid boundary to eliminate the need for 
a two-fluid model. Evaporation/condensa- 
tion at the liquid free surface was neglected. 
Since detailed measurements of the pressur- 
ant gas flow rate were not obtained during 
preliminary analyses that the flow rate was a 
constant value computed by dividing the 
total mass added during the run by the pres- 
surization time. Heat transfer rates to the 
tank wall and to the liquid /slush were taken 
from SLURP code analyses 4 , 

(Q w = 5.02 Btu/sec,Q 1/a =0.46 Btu/sec). The 


diffuser at the pressurant gas inlet was 
modeled as a porous media. 

Since all of the planned analyses were 
defined by axisymmetric boundary condi- 
tions for axisymmetric geometries, the first 
goal was to establish the validity of replacing 
a true three-dimensional model, (radial, 
azimuthal, axial), with a pie-sector model, 
(radial, 1, axial), to reduce the computational 
demands for completing the study. The first 
configuration under examination was ana- 
lyzed using a pie-sector model (10x1x10) and 
a fully three-dimensional model (10x36x10). 
Comparing the predicted velocity and tem- 
perature fields at a slice through a diameter 
of the fully three-dimensional model to those 
predicted by the pie-sector model revealed 
nearly identical solutions. To establish 
confidence that the pie-sector mesh was 
sufficiently refined, a convergence study was 
performed using a 15x1x15 mesh and a 
10x1x20 mesh. Although the finer meshes 
revealed some details missing from the first 
mesh, the predicted fields for the finer 
meshes did not qualitatively differ from those 
predicted by the first mesh. Since the liquid 
free surface was modeled as a solid boundary 
with a specified heat flux, the improved 
temperature gradient near the free-surface 
predicted using the finer meshes is not a 
significant factor for these analyses, but 
could become significant for analyses In 
which heat and mass transfer at the free 
surface is computed. It was therefore 
concluded that a 10x1x10 pie-sector mesh 
would be used for all subsequent analyses in 
this study. 

The second step in the evaluation 
process was to compare computational pre- 
dictions to experimentally obtained data. 

Since the pressurant Inlet pressure was not 
measured during the experiment, a value of 
50 psl was assumed for the Inlet pressure 
based on other known parameters In the 
experiment. The pressurant inlet mass flow 
rate is a key parameter which must be speci- 
fied to the computational model but which 
was not recorded during the experiment. The 
S-curvc of the pressure history displayed In 
Fig. 2 suggests that, contrary to Initial 
assumptions, the flow rate during the exper- 
iment was a function of time. A sequence of 
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analyses were performed during which the 
flow rate history was iteratively adjusted 
until the pressure history predicted by 
FLOW- 3D matched the history displayed In 
Fig 2. 

The difference between measured and 
computationally predicted total mass added 
was less than 15%. Although this difference 
Is larger than had been hoped for, several 
factors contributed to this difference. As 
detailed above. Important Input parameters, 
such as pressurant inlet pressure and flow 
rate history, were estimated or back- 
calculated because they had not been 
measured during the experimental effort. 
Given the level of uncertainty In the problem 
definition, the agreement between experi- 
mental and computational results was 
considered adequate to proceed with the 
Investigation. To eliminate this uncertainty 
In future efforts, a list of recommended 
measurements for future experimental 
studies Is provided In the concluding 
remarks of this paper. 


Parametric Study or Pressurization 
Process 

Table 1 presents a summary of the 
Input parameters and the performance pre- 
dictions for the twelve cases studied during 
this effort. An Initial pressure of 17.4 psi 
Indicates pressurization of liquid hydrogen 
whereas an initial pressure of 1. 1 psia Indi- 
cates that slush hydrogen was pressurized. 
For all analyses, the pressurant gas inlet 
pressure was specified to be 50 psi, the inlet 
temperature to be 307 °R, and the heat flux 
rates were specified to be Q w =5.02 Btu/sec 
and Qj/ S =0.46 Btu/sec. The pressurant 
mass flow rate for each case was assigned a 
constant value which is listed In Table 1. 

The analysis matrix defined In Table 1 repre- 
sents studies of different parametric effects 
which are described in detail in the balance 
of this section. 

The first study focused on the 
influence of initial temperature distribution 
within the ullage on the pressurization 
process. The ullage temperature distribution 
had not been measured during the experi- 
mental investigation so two bounding cases 


Table 1 Analysis Matrix (55.6% ullage) 



•Initial 

•Initial 

•Inlet 

••Press. 

••Mass 

••Flow 

••Final 

Case # 

Pressure 

Temp. 

Velocity 

Time 

Added 

Rate 

Pressure 


(psia) 

(°R) 

(ft/s) 

(s) 

(KW 

db m /s) 

(psia) 

1 

17.4 

60-40 

1.393 

33.0 

0.771 

0.0234 

50.2 

2 

1.1 

60-40 

1.393 

54.0 

1.262 

0.0234 

49.8 

3 

17.4 

250-40 

1.655 

26.0 

0.720 

0.0277 

49.8 

4 

17.4 

60-40 

1.655 

26.0 

0.720 

0.0277 

49.8 

5 

1.1 

60-40 

1.655 

43.0 

1.193 

0.0277 

50.4 

6 

1.1 

250-40 

1.655 

43.0 

1.193 

0.0277 

50.4 

7 

17.4 

60-40 

2.500 

15.5 

0.649 

0.0419 

49.5 

8 

1.1 

60-40 

2.500 

25.5 


0.0419 

50.1 

9 

17.4 

60-40 

3.500 

10.5 

0.612 

0.0583 

49.3 

10 

1.1 

60-40 

3.500 

17.0 

0.992 

0.0583 

49.2 

11 

17.4 

60-40 

5.000 

7.2 

0.600 

0.0832 

49.7 

12 

1.1 

60-40 

5.000 

11.5 

0.958 

0.0832 

49.3 


* specified value 

** computationally predicted value 
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were selected: the 60-40 °R distribution was 
selected to represent a fairly well mixed 
ullage whereas the 250-40 °R distribution 
was selected to represent a thermally strati- 
fied ullage. The temperature fields predicted 
at the end of the pressurization process for 
both temperature distributions and for both 
propellants (cases 3, 4, 5, 6) are displayed In 
Fig. 3. Review of Table 1 reveals that, given 
the same heat flux rate at the tank wall, the 
pressurization time appears to be Indepen- 
dent of Initial ullage temperature distribu- 
tion. Although the global pressurization 
performance as measured by total mass 
added and pressurization time Is unchanged 
by the different initial ullage temperature 
distributions, the detailed temperature pro- 
files do show some differences. In particular, 
the clustering of temperature contours at the 
surface of the slush hydrogen (for the low 
intlal pressure) Indicates that there may be a 
significant difference In heat transfer and 
mass transfer rates at the slush/vapor 
Interface than at the liquid/ vapor Interface. 
Since the current model uses a specified heat 
transfer rate at the free surface, and assumes 
zero mass transfer at the free surface, this 
hypothesis Is unprovable with the present 
modeling capability and remains a topic for 
future study. 

The relationship between pressuriza- 
tion time and pressurant flow rate revealed 
by the matrix of cases In Table 1 is displayed 
In Fig. 4. Although the decrease in pressur- 
ization time with Increasing mass flow rate Is 
expected, it should be noted that the rela- 
tionship is not a simple Inverse proportion- 
ality. Figures 5 and 6 present one of the 
most interesting results of the study. The 
computationally predicted relationship 
between pressurization time and the total 
pressurant mass addition required to 
pressurize the tank to 50 psi Is displayed in 
Fig. 5. The relationship between pressurant 
mass flow rate and the total mass addition 
required Is displayed In Figure 6. The two 
relationships are consistent and lead to the 
conclusion that a more rapid pressurization 
requires less pressurant mass than a slow 
pressurization. Although this conclusion 
may be counterintuitive, the explanation lies 
In the Imposed thermal boundary conditions. 
A longer pressurization time, coupled with 


the constant heat flux specified at the tank 
walls and at the liquid/slush Interface, 
results In a greater cooling of the ullage. At 
the lower vapor temperature, a greater vapor 
mass Is required to reach the desired 
pressure level. Review of the SLURP study 
report 4 reveals data to support this 
conclusion and a comment on this trend by 
those researchers. Therefore, the ability of 
the FLOW-3D model to correctly predict this 
trend lends additional confidence to the 
validity of the computational model, and 
brings additional attention to this interesting 
and potentially important phenomena. 

The dependence of velocity and 
temperature fields In the ullage at the end of 
pressurization on the pressurant mass flow 
rate was examined by extracting appropriate 
data from the results generated by the anal- 
ysis matrix defined In Table 1 . All the 
analyses examined for this study began with 
the 60-40 °R ullage temperature distribution 
and proceeded from an initial saturation 
pressure to a final pressure of 50 psi. Fig- 
ure 7 presents computed velocity fields in the 
ullage for pressurization of liquid hydrogen 
and Fig. 8 presents the corresponding 
temperature fields. For the range of Inlet 
mass flow rates examined, the entering gas 
Jet never fully penetrated into the ullage 
volume. As the flow rate Increases a vortex is 
established at the top of the tank near the 
inlet diffuser. The temperature near the free 
surface varies from 50.5 °R to 57.4 °R for the 
range of pressurant mass flow rates studied. 
Although the presence of the jet Is apparent 
In the temperature contours, the bulk of the 
ullage space is well represented by an axial 
temperature distribution with little radial 
dependence. 

The velocity fields computed for pres- 
surization of slush hydrogen are presented in 
Fig. 9 and the corresponding temperature 
fields In Fig. 10. Comparing Fig. 9 to Fig. 7 
reveals that substantially different flow pat- 
terns are established during slush pressur- 
ization than during liquid pressurization. At 
the lower flow rates, the incoming Jet again 
fails to fully penetrate the ullage and estab- 
lishes a vortical flow near the Inlet. At the 
higher flow rates, the Incoming pressurant 
Jet Is seen to fully penetrate the ullage 
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volume and to establish a variety of flow 
patterns. At an Intermediate flow rate, the Jet 
penetrates all the way to the slush surface 
and then follows a serpentine path establish- 
ing two side-by-side vortices. At the highest 
flow rate, the Jet penetrates to the slush 
surface and moves radially outward along the 
surface toward the tank wall. Although two 
vortices are again formed, they are arranged 
one-above-the-other Instead of side-by-side. 
Examining Fig. 10 shows that the variety of 
flow patterns predicted for slush pressuriza- 
tion is reflected in the variety of ullage 
temperature distributions. At the higher flow 
rates, the uniformity of temperature with 
radial location observed for liquid pressur- 
ization is replaced with highly radially 
dependent temperature distributions. The 
predicted ullage temperature near the slush 
surface increases from 169 °R for the lowest 
Inlet mass flow rate to 286 °R for the highest 
flow rate, as compared to approximately 
60 °R for the liquid cases. 

Comparing the predicted temperature 
and velocity fields for pressurization of liquid 
hydrogen to those predicted for slush hydro- 
gen revealed some significant differences. 
These differences in turn raise concern about 
some of the modeling assumptions used to 
perform the analyses. The highly stratified 
ullage temperature distributions produced by 
pressurization of liquid hydrogen lead to 
relatively low gas temperatures near the 
liquid/vapor interface which, when coupled 
with the nearly stagnant state of gas near the 
free surface, seem to support the modeling 
assumption of constant heat flux with 
negligible mass transport at the liquid/vapor 
interface. In contrast, high gas temperatures 
predicted during pressurization of the slush, 
coupled with significant flow velocities 
predicted near the free surface, lead one to 
seriously question this modeling assump- 
tion. Since evaporation or condensation can 
have a significant effect on the tank Internal 
pressure, caution should be used In accept- 
ing computational predictions which do not 
include these effects. To definitively resolve 
this concern requires modeling capabilities 
beyond the those of the present 
computational tool. 


The evolution of temperature and 
velocity fields observed during the Inlet mass 
flow rate dependence study naturally lead to 
a curiosity about the dependence of the tank 
pressure history on inlet mass flow rate. 
Figure 1 1 displays the tank pressure histo- 
ries predicted for pressurization of liquid 
hydrogen. The tank pressure histories 
predicted for pressurization of slush are 
displayed In Fig. 12. Although the trends are 
not surprising, it should be noted that the 
rate-of-change of pressure Is not constant 
during the pressurization process but rather 
increases slightly with increasing elapsed 
time. A desire to understand this behavior 
lead to the modeling and dimensional 
analysis described in the next section. 


A Simplified Model and Non-Dimensional 
Characterization 

If multidimensional effects of the 
pressurization process are neglected, and the 
ullage is modeled as a lumped thermody- 
namic entity, a simplified analytical model of 
the pressurization process can be readily 
developed. Conservation of thermal energy 
can be enforced by satisfying a balance 
equation written for the ullage volume: 

w d(pT) 

d ~ = m g h g - Qw- Qys m 

where Q w includes any discrete "heat leaks” 
Into the ullage as well as convective heat 
transfer from the ullage to the tank wall. The 
product of density with temperature in the 
derivative can be replaced by an expression 
including pressure by assuming that the 
ullage can be modeled as an ideal gas 
(pT=p/R). 

dt = ^(“gV Qw Ql/s) 121 

This form clearly shows that, for constant 
specified heat flux at the boundaries and 
constant pressurant mass flow rate, this 
lumped model predicts a constant rate-of 
change of tank internal pressure with time. 
This is in qualitative agreement with the 
near-constant pressurization rate predicted 
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by the computational model. The difference 
Is due to the multidimensional effects Incor- 
porated Into the computational model which 
cannot be Included in the simple lumped 
model. 

Since one goal of the modeling and 
analysis Is to develop correlating parameters 
for the pressurization process. It is important 
to select appropriate normalizing or non- 
dlmenslonalizlng scales. One logical time 
scale is the elapsed time required to complete 
the pressurization process. 


Rearranging, and substituting appropriate 
values into this equation, produces the 
following result. 


^2 Pf* - pr 
IPf-PiJ dt = tf'-tj* 


1 -0 
1 -0 


= 1 


171 


Therefore, for any pressurization analysis 
with specified constant mass inflow rate and 
constant heat flux boundary conditions: 



[ 8 ] 



The total change in pressure to be accom- 
plished by the pressurization process was 
selected to normalize the pressure rise from 
pressurization initiation. 


Since this result is independent of specific 
values for any of the parameters required for 
the computation, the goal of a "universal" 
relationship for all pressurization processes 
has been satisfied (for the lumped parameter 
model). 


P - Pi 
Pf-Pi 


[41 


It was hoped that this combination of charac- 
teristic scales would produce a nondimen- 
sional representation of the pressurization 
process which would produce a single 
relationship for all the pressurization cases. 
Combining the nondlmenslonal variables 
defined in Eqs. 3 and 4 produces the 
following expression for the rate-of-change of 
nondimensional pressure with normalized 
time. 


d 

dp* 

P * Pi I 

LPf-Pi J 

dt* “ 

d 

t 1 
. fc f. 



Pf-Pi 


d£ 

dt 


[5] 


Since the computational model 
includes multidimensional effects which 
cannot be included in the lumped parameter 
model, the validity of dp*/ dt* as a correlating 
parameter for pressurization analyses must 
be evaluated against the results from the 
analysis matrix defined in Table 1. The tank 
pressure history data previously displayed in 
Figs. 1 1 and 12 were reorganized using the 
new dimensionless variables to produce 
Fig. 13. Careful examination reveals that, 
although dp*/dt* is not precisely equal to a 
constant value of one, reorganizing the data 
using p* and t* does collapse all the data into 
a single curve. The deviation of this curve 
from a straight line is similar to that seen in 
the dimensional plots and is due to multidi- 
mensional effects. 


Recalling that dp/dt Is constant for constant 
heat flux and flow rate boundary conditions, 
integrating this expression over the elapsed 
time required for pressurization yields the 
following result. 


Pf* - Pi* = 


*f dp 

.Pf-Pi dt. 


(tf* - V) 


16) 


Summary and Concluding Remarks 

The proposal to use slush hydrogen as 
the propellant storage medium for the NASP 
presents a formidable challenge to the design 
and efficient operation of the propellant 
management system. Previous studies of 
propellant tank pressurization prior to main 
engine ignition have primarily involved 
experimental measurements In model 
tanks 1,2 or analytical and computational 
models based on lumped or one-dimensional 
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models 4 . The present study was initiated: to 
develop a capability for multidimensional 
modeling of the propellant tank pressuriza- 
tion process prior to main engine ignition; to 
evaluate the performance of, and identify 
limitations in, the model and; to begin eval- 
uation of the parametric dependence of the 
pressurization process on the system initial 
conditions and boundary conditions. 

After a review of available computa- 
tional tools, the FLOW-3D code was selected 
as the most promising candidate for model- 
ing pressurization and expulsion of slush 
hydrogen from storage tanks. Code perfor- 
mance was evaluated by comparing compu- 
tational predictions to measurements 
obtained during the pressurization of an 
experimental simulation of a NASP propellant 
tank 1 . The analyses demonstrated that, 
despite significant uncertainties in problem 
definition due to incomplete specification of 
the experiment conditions, a pie-sector 
(axisymmetric) computational model was 
able to replicate the tank pressure history 
and to predict the total pressurant mass 
required for the pressurization to within 15% 
of the measured value. Based on this 
evaluation, it was concluded that FLOW-3D 
is a suitable tool from multidimensional 
modeling of the pressurization process. 

Typical of many efforts to verify 
computational tools by comparison to exper- 
iment data not specifically obtained for this 
purpose, the performance evaluation was 
complicated by the absence from the experi- 
ment report of key parameters required to 
define the model. To help e liminat e 
uncertainties in model definition, future 
experiments should include measurement of: 
the history of pressurant inlet pressure, 
temperature, and mass flow rate; the initial 
temperature distribution within the ullage; 
the history of heat loss to the tank wall and 
to the slush/vapor interface and; the rate of 
evaporation or condensation at the free 
surface. Measurements necessary to 
evaluate code performance should include 
the pressure history in the tank and the 
transient velocity and temperature fields in 
the ullage during pressurization. Although it 
is recognized that the distinction between 
model definition and performance evaluation 


measurements is tenuous, and that some of 
the requested measurements are difficult to 
obtain, the lists have been assembled as a 
starting point for consideration by future 
efforts. 

The dependence of the pressurization 
process on initial ullage temperature 
distribution was examined by comparing 
predictions for a 60-40 °R distribution, 
representing a well mixed initial condition, to 
those for a 250-40 °R distribution, represent- 
ing a thermally stratified initial condition. 

For the pressurant mass flow rates consid- 
ered and the heat flux rate specified at the 
ullage boundary, the global performance as 
measured by elapsed time required for the 
pressurization and the total pressurant mass 
added appears to be nearly Independent of 
the initial ullage temperature distribution. 
Review of the predicted velocity and 
temperature fields near the free surface at 
the end of pressurization indicates that the 
assumption of no mass transfer at this Inter- 
face merits future scrutiny. Further, the 
specification of a single heat flux rate at the 
tank wall for all cases is also suspect and 
should be replaced by a more accurate model 
when possible. 

A study of the effect of pressurant 
mass flow rate on the pressurization process 
revealed an interesting trend as well as 
potentially significant differences between 
pressurization of liquid hydrogen and 
pressurization of slush hydrogen. For both 
liquid and slush cases, the pressurant mass 
required to reach the desired pressure level 
decreases with increasing pressurant mass 
flow rate. This Is in accord with experiment 
data 1 and is discussed in detail in the 
preceding section. Review of the temperature 
and velocity fields predicted for the slush and 
liquid cases during this study revealed 
interesting differences between the cases. At 
comparable pressurant mass flow rates, the 
inlet Jet fully penetrates the ullage of the 
slush cases whereas it does not for the liquid 
cases. This difference leads to significantly 
different flow and temperature fields in the 
ullage. The temperature field predicted for 
the liquid cases has little radial dependence 
and the one-dimensional model provided by 
the SLURP code 4 is probably reasonable. In 
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contrast, the temperature fields predicted for 
pressurization of slush have a very signifi- 
cant radial dependency and will not be well 
represented by a one-dimensional model. 
Further, at the higher mass flow rates, the 
temperatures and velocities near the free 
surface of the slush are significantly higher 
for the slush cases than for the liquid cases. 
The low values predicted for the liquid cases 
suggest that the assumption of no mass 
transfer at the free surface may be reason- 
able whereas the high values for the slush 
case cast serious doubt on the validity of 
such an assumption. 

A simplified model of the pressuriza- 
tion process was developed In search of a 
nondlmenslonal characterization of the 
process. Based on this model, a dimension- 
less time was defined as elapsed time from 
pressurization Initiation normalized by the 
time required to complete pressurization. A 
dimensionless pressure was defined by the 
increase In pressure from Initial pressure 
divided by the total pressure rise to be 
accomplished. Using these dimensionless 
parameters, it was demonstrated that the 
tank pressure history for all of the cases In 
this study collapse Into a single curve. 
Although the relationship revealed by this 
study lacks predictive capability since the 
dimensionless parameters contain variables 
which are not known a priori, it Is valuable 
for unifying the presentation of data from 
experimental and computational studies of 
the pressurization process. 
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Figure 1 Schematic of experiment tank 


Tank pressure (psi) 



Figure 2 Pressurization of tank containing 
liquid hydrogen (55.6 % ullage) 
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(c) Case # 5 (d) Case # 6 


Figure 3 Temperature distribution, at the end of pressurization (m — 0.277 lbm/s) 
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Figure 4 Pressurization time .vs. 
Mass flow rate 


Mass added (Ibm) 



Figure 6 Mass 


Mass added (Ibm) 



Figure 5 Mass added .vs. Pressurization time 



led .vs. Mass flow rate 


2 










(a) Case # 1 


(b) Case # 7 




Figure 7 Velocity fields for pressurizing liquid hydrogen 
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(a) Case # 1 


(b) Case # 7 



(c) Case # 9 


(d) Case # 11 


Figure 8 Temperature fields for pressurizing liquid hydrogen 


14 








(a) Case # 2 


(b) Case # 8 



(c) Case # 10 (d) Case # 12 


Figure 9 Velocity fields for pressurizing slush hydrogen 
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(c) Case #10 


(d) Case #12 


Figure 10 Temperature fields for pressurizing liquid hydrogen 
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Figure 1 1 Tank pressure history 
(liquid hydrogen) 
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Figure 12 Tank pressure history 
(slush hydrogen) 



Figure 13 Dimensionless pressure .vs. 
normalized time 
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Abstract 

A multi-dimensional computational model of the 
pressurization process in liquid/slush hydrogen ^ is 
developed and used to study the influence of heat flux rates 
at the ullage boundaries on the process. The new model 
computes these rates and performs an energy hainner for 
the tank wall whereas previous multi-dimensional mod e ls 
required & priori specification of the boundary heat flux 
rates. Analyses of both liquid hydrogen and slush 
hydrogen pressurization were performed to expose 
differences between the two processes. Graphical displays 
sxc presented to establish the dependence of pressurization 
time, press urant mass required, and other parameters of 
interest on ullage boundary heat flux rates and press urant 
mass flow rate. Detailed velocity fields and temperature 
distributions are presented for selected cases to further 
illuminate the details of the pressurization process. It is 
demonstrated that ullage boundary heat flux rates do 
significantly effect the pressurization process and tha t 
minimizing heat loss from the ullage and maximizing 
pressurant flow rate minimizes the mass of press urant gas 
required to pressurize the tank. It is further demonstrated 
that proper dimensionless scaling of pressure and time 
permit all the pressure histories examined during this study 
to be displayed as a single curve. 


Nomenclature 

A area 

h heat transfer coefficient 

MC mass-specific heat product 

P pressure 

Pi initial pressure 

Pf final pressure 

P dimensionless pressure 

t time 

tf pressurization time 

t* dimensionless time 

T temperature 

T w tank wall temperature 

Q w rate of heat transfer from ullage to tank wall 

Qi rate of heat transfer from ullage to liquid/slush 
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Introduction 


we past several years, the pressurization process 
in cryogenic propellant tanks has been the focus of several 
analytical 1 ’ 2 , experimental 2 - 3 , and numerical studies 4 . 
Design and efficient operation of liquid/slush hydrogen 

propellant tanks (such as proposed for the National 

Aerospace Plane), requires the ability to accurately model 
this process. The use of slush hydrogen propellant makes 

this a particularly challenging task because it requires the 

modeling of transient flmri m Mttl « • 


processes in the presence of all three phases of matter. The 

flow field and the temperature field in a tank daring 
pressurization are multi-dimensional and it is likely that the 
tank will be subject to multi-dimensional environmental 
infl u en ces. Existing models are fiiMamwitcn y ^ 
dimensional and based on simplifying assumptions and 
experimentally derived correlations. Althou gh one- 
dimensional models may be satisfactory for p reliminar y 

design studies, it is not clear that they are adequate as 

design tools for future vehicle development Sasmal, 
Hochstein, and Wendl 4 reported a multi-dimensional 
model for such a process using the FLOW-3D code 5 and 
demonstrated its accuracy by comparing computational 
predictions to experimentally ob tained data They 
simplified the heat transfer models at the ullage boundaries 
by specifying constant heat flux rates predicted by a one- 
dimensional model, (the SLURP code) 1 , to study the 
psraroetnc dependencies of the pressurization process. 
Although a reasonable first approximation, this assumption 
is suspect because the process itself is transient and is 
influenced by several parameters which have values that 
change during the pressurization process. 

The present study was undertaken to more accu ratel y 
model the heat transfer at the ullage boundaries and to 
investigate its effect on the pressurization process. 

Although the predictions of the present model are in 
qualitative agreement with previously published data 4 , the 
quantitative differences produced by the improved transient 
heat flux rates at the ullage boundaries clearly indicate the 
strong dependency of the pressurization process in liquid 
and slush propellant tanks on these quantities. 
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Propellant Tank Configuration 

The geometric configuration used in the present study 
is a partially filled 5 ft diameter spherical aluminum tank 
with a wall thickness of 03 inch. A schematic of this 
configuration is presented in Fig. 1. 



llie tank is pressurized by injecting gaseous hydrogen 
through a 1 ft diameter diffuser located at the top of the 
tank until the tank pressure reaches 50 psi. The initial 
vapor pressure is 17.4 psi when pressurizing liquid 
hydrogen and 1.1 psi when pressurizing slush hydrogen. 
These specifications were adopted from an experimental 
study 0 so that computational predictions could be 
compared to experimentally obtained data for the 
pressurization process. 


Computational Model 

The FLOW-3D code 5 was selected as the 
computational tool to provide multi-dimensional modeling 
of the pressurization process. Since the configuration of 
interest is an axisymmetric geometry, and the present study 
is focused on axisymmetric boundary conditions, a pie 
sector model (10x1x10 mesh) was selected to reduce the 
computational demands far completing the study. Pie- 
sector model performance, computational grid 
convergence, and the suitability of FLOW-3D for modeling 
the pressurization process in liquid/slush hydrogen 
propellant tanks has already been documented*. A 
previous pressurization study 4 enforced constant heat flux 
rates at the ullage boundaries that were obtained from a 
one-dimensional model 1 . More accurate modeling of these 
boundary conditions was achieved in the present study by 
dividing the tank wall into azimuthal strips and identifying 
each as a separate obstacle. The code provides an option to 
model each obstacle as a thermal lump and can compute 
the temperature of each obstacle using the expression: 

Mc ^=hA( T -T.) („ 


where MC is the mass-specific heat product, A is the 
obstacle surface area, h is the heat transfer coefficient, T 
is the obstacle temperature, and T is the fluid temperature 
in the adjacent cell. The code computes values for h based 
on correlations for natural convection, forced convection 

and conduction within the fluid and then uses the largest ’ 

value when computing the wall temperature 


Analyses, Results, and Discussion 


The primary goal of the present study is to examine 

the dependence of the pressurization process on the heat 

transfer rates from the ullage to the tank wall and to the 

liquid/slush pooL Hie traditional approach to 
multivariable problems, (varying one parameter at a time) 

was used to study the sensitivity of the pressurization 

process to "assumptions" as well as traditional "control" 
variables. Where available, the initial and boundary 
conditions specified for the analyses were those measured 
dunng the corresponding experimental study 8 . Although a 
detailed axial temperature distribution was not measured 
<hmng the experimental study, the initial temperature in the 
image is approximated by a linear variation from 250° R at 
the inlet diffuser to 40° R at the phase interface . The 
initial tank wall temperature was therefore specified to 

vary linearly from 220° Rat the diffuser to 40° Rat the 

interface. The inlet gas temperature was a constant 307° R 
and tire inlet pressure was maintained at 50 psi. A constant 
rate of heat flux per unit area from the environment to the 
tank (0.00884 Btu/sec-ft 2 ) was specified. To simplify the 
present study, no evaporation/condensation at tire 
liquid/slush-vapor interface was assumed. All analyses 
were performed for an ullage volume fraction of 0356. 
Using this model, analyses were performed for 
pressurization of a liquid and a slush propellant tanir to 
reveal differences in their response to similar boundary 
conditions. 


To expose the sensitivity of the pressurization process 

to tire heat transferrate from the ullage to tire tank wall and 
to the liquid/slush-vapor interface, a sequence of analyses 
were performed. The first analysis used tire same heat 
transfer rates computed by a one-dimensional SLURP 
model 1 that were used in the previous study 4 . The second 
analysis retained the same liquid/slush-vapar interface heat 
transfer rate, but enforced an adiahatic boundary at the tank 
wall. Although this is not an improvement in the fidelity 
of the model, it is useful as an extreme condition for 
examining the sensitivity of tire pressurization process to 
this parameter. The third analysis used the obstacle heat 
storage model for the tank wall and a specified temperature 
boundary condition to model the liquid/slush-vapor 
interface. The linear temperature variation 

above was used to initialize the tank wall obstacles and the 

bqui d/slush-vapor interface was specified to be a constant 
40“ R. Figure 2 displays the relationship between required 
pressurant mass and pressurant mass flow rate for 




Figure 2 Mass added vs. Mass flow rate (liquid) 
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Figure 3 Mass added vs. Mass flow rate (slush) 


pressurization of a liquid propellant tank using the three 
different heat transfer boundary conditions. Figure 3 
displays the corresponding information for pressurization 
of a slush propellant tank. The cases with heat transfer rate 
computed by FLOW-3D have higher heat losses than the 
other cases and reveal a decrease in required press urant 

mass with decrease in the rate of heat loss from the ullage. 

With no heat loss to the tank wall, and 0.46 Btu/sec heat 
loss rate to the liquid/slush hydrogen, press urant mass 
required was neatly constant for all mass flow rates. 

Results from analyses with heat transfer to the tank wall 
show a decreased pressurant mass requirement with 
increasing mass flow rate. Comparing cases with si milar 
boundary models, the slush hydrogen consistently requires 
more pressurant mass than the liquid to reach 50 psi. 
Review of the results from these analyses leads to the 
conclusion that the more detailed model consistently 


predicts a larger required pressurant mass and is therefore 
conservative for design purposes. 

The tank wall temperature distribution predicted at 
the end of pressurization for different inlet mass flow rates, 
using the obstacle heat storage model, are presented in 
Figs. 4 and 5 for liquid and slush hydrogen respectively. 

Ditanc* from booam ot •» urk (IQ 



Figure 4 Wall temperature distribution (liquid) 



Figure 5 Wall temperature distribution (slush) 


The final tank wall temperature distributions show that 
higher temperatures are predicted for the slush case than 
for the liquid case. This difference is believed to be due to 
the longer pressurization time required for the «!»■<* case 
and a significantly different flow pattern in the »tia g » for 
the two cases. Figures 6 and 7 display the predicted 
temperature distributions along the tank centerline at the 
end of the pressurization period for liquid and for «in«h 
hydrogen respectively. These figures clearly display the 
extent of penetration of the hot gas into the cold ullag e and 
reveal significantly greater penetration for the slush case 
than far the liquid case. This interpretation is further 
substantiated by detailed velocity and temperature fields 
discussed later. 

The influence of pressurant mass flow rate on the rate 
of heat loss from the ullage to the tank wall is displayed in 
Fig. 8 for pressurization of liquid hydrogen and Fig. 9 for 
pressurization of slush hydrogen. These profiles were 
predicted using the obstacle heat storage model. For both 
cases, the rate of heat loss increases with increasing time 
and increases mare rapidly at higher inflow rates. This is 
reasonable, since higher mass flow rates should result in a 
higher heat transfer coefficient at the wall. 
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Figure 6 Ullage temperature distribution (liquid) 



Figure 8 Heat transfer rate from ullage to tank wall (liquid) 



Figures 10 and 1 1 (liquid and slush) display a decrease in 
pressurization time with increase in mass flow rate as one 
would expect from a cursory analysis. The increased heat 
loss rates computed using the obstacle heat storage model 
as compared to the other models result in increased 
pressurization times. 



Figure 7 Ullage temperature distribution (slush) 



Velocity and temperature fields predicted for slush 
and for liquid hydrogen pressurization using the obstacle 
heat storage model are depicted in Figs. 12-15. These 
fields are different from each other and from thn^f 
predicted using SLURP-based constant boundary heat flux 

rates 4 . Comparing the predicted fields reveals significant 


4 










figure 12 End of pressurization temperature fields (liquid) Figure 13 End of pressurization velocity fields (liquid) 



figure 14 End of pressurization temperature fields (slush) Figure 15 End of pressurization velocity fields (slush) 
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differences due to changes in heat flux rates and depending 
on whether slush or liquid hydrogen is being pressurized. 
As already noted, the hot pressurant gas penetrates the 
ullage more deeply when pressurizing slush than when 
pressurizing liquid hydrogen. 

The final task undertaken during this study was to 
examine the predicted pressure history using the 

dimensionless pressure, P* = IP ~ Pl * and 

(P» - Pi) 

dimensionless time, T* = — , developed during the 

previous study 4 . Figure 16 displays the pressure history 
predicted using the SLURP-generated constant boundary 
heat flux rates. Fig. 17 displays the history predicted using 
the adiabatic tank wall boundary, and Fig. 18 displays the 
history predicted using the obstacle heat storage model. 
When overlaid, these curves are close enough to be 
considered the same. Therefore, the pressure history when 
plotted using the dimensionless parameters is not 
significantly influenced by any of the parameters e xamine 
during the present study. 

Summary and Conclusions 

A multi-dimensional computational model of the 
pressurization process in liquid/slush hydrogen tanks was 
developed using the FLOW-3D computer code. In contrast 
with previously published multi-dimensional models, the 
heat .flux rates from the ullage to the tank wall and to the 
liquid/slush are computed by the model rather than 
requiring & priori specification. The tank wall is modeled 
using a lumped heat storage model applied to azimuthal 
strips of the wall and the heat flux rates are computed from 
standard correlations. 

The new model was used to study the influence of 
ullage boundary heat flux rates on the pressurization 
process. Analyses of both liquid hydrogen and slush 
hydrogen pressurization were performed to expose 
differences between the two processes. Predictions from 
the new model were compared to previously published 
predictions which used constant heat flux rates predicted 
by a one-dimensional model and to analyses based on an 
adiabatic interface between the ullage and the tank wall 
For each of these models, analyses were performed over a 
range of pressurant mass flow rates. The pressurization 
predictions from these analyses are summarized in 
graphical displays to establish the functional 
interdependence between the parameters. 

It is demonstrated that the ullage boundary heat flux 
rates do significantly effect the pressurization process. The 
pressurant mass requirement is shown to increase with 
increasing heat transfer rates at the ullage boundaries and 
to also increase with decreasing mass flow rate. It is 
shown that the entering pressurant gas more deeply 
penetrates the ullage when pressurizing slush than when 
pressurizing liquid and 
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that this results in a significant change in the temperature 
distribution in the ullage and the tank wall. The heat 
transfer rate from the ullage to the tank wall is shown to 
increase during the pressurization process and with 
increasing mass flow rate. Detailed velocity fields and 
temperature distributions are presented for selected cases to 
further illuminate the details of the pressurization process. 
Previously defined dimensionless pressure and time 
parameters are used to organize all of the pressure history 
profiles generated during this study. It is demonstrated that 
proper dimensionless scaling of pressure and time permit 
all the pressure histories examined to be displayed as a 
single curve. 
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